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ABSTRACT 

Context. Studying circumbinary disks is critical to understanding the formation mechanisms of binary stars. While optical or mid-infrared 
images reveal the scattered emission, millimeter observations provide direct measurements of the dust thermal emission. 
Aims. We study the properties of the circumbinary disk around the well-known, multiple young stellar object GGTau with the highest possible 
sensitivity and spatial resolution. 

Methods. We mapped the continuum emission of GGTau at 267 GHz using the IRAM Plateau de Bure interferometer equipped with upgraded 
receivers and LO systems. An angular resolution of 0.45" x 0.25" was achieved, corresponding to a linear resolution of 65 x 35 AU. 
Results. The GGTau A circumbinary disk is observed as an extremely clearly defined narrow ring. The width of the ring is not resolved. 
Emission from the central binary is detected and clearly separated from the ring: it coincides with the GG Tau Aa position and may therefore 
trace a circumstellar disk around this star. The mass ratio of the circumbinary to circumprimary material is ~80. 
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1. Introduction 

GG Tau is a well-known quadruple system located in the Taurus 
star-forming region. It is composed of two binaries, the north- 
em one (GG Tau A) having been the subject of considerable 
attention during the past decade. This system, classified as a 
classical TTauri star, is a close binary with a separation of 
~35AU (for an assume d distance of 140 pc). Its age is esti- 
mated to be 1.5 Myrs dWhite et all 1 1999b . GGTau A is sur- 



rounded by a circumbinary disk, which has been extensively 
studied by means of its visible or mid-infrared sca t erred emis- 
sion (e.g. Roddieretal. 1996; Silber et al.' '200d iKrist et al 



l2002l;ElcCabe et al.. 2002: Krist et al. .2005.) and its direct ther- 



mal dust emission in the mill imeter domain (e.g. lDutrev et al 
1994 : Guilloteau et al. 1999 ^. The higher angular resolution 
data clearly illustrate that the inner ~ 150 AU region of the disk 
has been emptied by tidal interactions with the two stars, mak- 
ing GG Tau one of the text-book example of circumbinary disk 
in young stellar syst e ms. 

h999h information 



Guilloteau et al. 



from CO and HCO lines, demonstrating that the disk is in 
Keplerian rotation (the redshifted emission being on the west- 
em side) around a ~ 1 .3 Mq system. These IRAM Plateau de 
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* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



Bure observations also indicated that the circumbinary disk is 
composed of a ring with sharp edges (~70% of the mass), sur- 
rounded by a more extended, c older disk. The to t al ring n-disk 
mass is ~O.12M0. Interestinglv. lGuilloteau et al.l (Il999h noted 



the continuum 1 .4 mm emission also revealed the presence of 
material within the ring, suggesting that there are circumstellar 
disks around the individual sources. 

In this paper, we present new observations of the millimeter 
continuum emission of the GG Tau circumbinary disk, which 
provides a significant improvement in angular resolution and 
sensitivity compared to the Guilloteau et al. (1999) data and 
places new constraints on the properties of the circumbinary 
and primary disks. 

2. Observations 

During the winter 2007/2008, new local oscillator components 
for the 1 mm receiver band were installed on the Plateau de 
Bure Interferometer, enabling observations up to a frequency 
of 270 GHz (instead of 250 GHz previously). As part of the 
science verification program of this upgrade, GG Tau was ob- 
served in January-February 2008. The receivers were tuned in 
single side band at 267.2 GHz. Only one polarization chan- 
nel was used, with a correlator setup covering a bandwidth of 
2 GHz to enable the simultaneous observation of the HCN(3- 
2) and HCO^(3-2) lines (which data will not be discussed in 
this article). 
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Fig. 1. Real part of the visibility vs baseline lengths. The points and error bars represent the measured values, and the curves are 
the best fits of a geometrical ring (plain curve; see Table [T) and of a ring plus a point source (dashed). The visibilities have been 
corrected from the disk orientation (PA = 7°) and inclination (/ = 37°) by compressing the V value by cos(/) prior to circular 
averaging. 



The source was observed with five antennas in the C config- 
uration (baseline length from 25 to 175 m) and with six anten- 
nas in the A configuration (baseline length from 130 to 760 m). 
Both observations were done in track-sharing mode with two 
other sources, resulting in a six-antenna-equivalent total on- 
source integration time of 4 hours. The atmospheric conditions 
were excellent, with a water vapour content of 0.5-1 mm. This 
translated into a ~200 K system temperature around transit. 

The data processing was performed with the GILDASQ 
software. Standard phase and amplitude calibration was ap- 
plied to the data, using regular observations of the B0548-(-134 
quasar. The atmospheric phase noise was 15^0° for the C 
configuration and 30-50° for the A configuration. The ab- 
solute flux calibration was performed using observations of 
MWC 349. Since this project was among the first Plateau de 
Bure observations at that high frequency, the flux density of 
MWC 349 was known only with a limited accuracy. We used an 
extrapolation of the measured flux densities at lower frequen- 
cies. We carefully checked the relative calibration of the two 
tracks by comparing the GGTau flux density on the overlap- 
ping baselines. We estimate the final accuracy of the absolute 
flux density to be ~ 20%. 

The data were imaged using robust weighting, resulting in 
one of the highest angular resolutions achieved with the Plateau 
de Bure Interferometer, 0.45x0.25" (at PA 19°). At the distance 
of 140 pc, this corresponds to a linear resolution of 65 x 35 AU. 

3. Results 

The image reveals a very clear ring structure as well as an un- 
resolved source located close to the ring center We performed 
a fit of a geometrical ring of uniform brightness to the data, 
whose results are given in Table [H They indicate a narrower 
ring than previous studies (IGuilloteau et al ] fl999h . because of 
the increase in angular resolution. Fig. [Tldisplays the real part 



Ring observational parameters 





This work 


Guilloteau et al. 1999 


Frequency (GHz) 


267 


220 


Ang. resolution(") 


0.45 X 0.25 


0.88x0.56 


Flux (mjy) 


760 + 10 


544 ±3 


Inner radius (") 


1.38 ±0.01 


1.27 ±0.01 


Outer radius (") 


1.84 + 0.01 


1.90 ±0.01 



Table 1. Parameters of a geometrical ring fitted to the uv data, 
and comparison with the results obtained with previous data. 



' See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about GILDAS. 



of the circularly averaged visibilities as a function of the base- 
line length together with the visibility curve of the fitted ring, 
showing the excellent agreement with the data. A marginally 
tighter fit can be obtained by fitting a ring and a point source 
(Fig-ID dashed curve). 

In addition to standard imaging, we used a more sophis- 
ticated technique to optimize the imaging process for this 
dataset. The dirty beam of the observations had somewhat high 
sidelobes (at about 25% level), owing to the poor uv coverage 
of the compact C configuration. As a consequence, the emis- 
sion within the (bright) ring was difficult to deconvolve prop- 
erly. To overcome these difficulties, we subtracted from the vis- 
ibilities the fitted geometrical ring (Table [TJ and imaged and 
deconvolved separately the ring and the residual data. We then 
added the two deconvolved maps to obtain the final image pre- 
sented in Fig. |2l The advantage of this method is the level of 
sidelobes caused by the bright ring structure is reduced and the 
dynamical range of the image is thus increased, especially in 
the inner part of the ring. On the final image, the peak bright- 
ness temperature is ~4 K, while the rms noise of the map is 
1.6 mJy/beam (0.24 K), i.e. the dynamical range is > 15. 
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Fig. 2. Left: 267 GHz continuum image of GGTau. The angular resolution is 0.45" x 0.25" at PA 19°, and the contour spacing 
is 3.1 mJy/beam (0.47 K, about 2 cr). The map is centered at a = 04 : 32 : 30.356, 5 = 17 : 31 : 40.48 (cross, J2000). The 
GG Tau B binary is located ~ 1 0" towards the south, outside the plotted area. Top right: zoom of the inner 1 " and position of the 
GG Tau A binary stars as of 2008. Bottom right: best physical ring model (see Table|2]and text). Angular resolution and contour 
spacing are as in the observations. 



4. Data analysis 

4.1. The ring 

To analyse the continuum emission of the GG Tau circumbi- 
nary ring, we fitted a physical model of a ring (in contrast to the 
geometrical model mentioned in the previous section). From 
a parametric density and temperature description, the model 
solves for the radiative tranfer to predict the observed emission, 
and performs a global minimization directly in the uv plane 
(see e.g. , Pietu et al. 2007, for a full description of the method). 
Here, the ring is assumed to be circular with a constant width 
and density. Its edges are modelled as Gaussians decreasing 
from the ring plateau. The model also includes a standard disk 
surrounding the ring, but its contribution is found to be neg 



ligible in reproducing these continuum data. ^ Guilloteau et al 



1999h showed that ~ 30% of the mass lies in the outer disk. 



but our data lack short baseline measurements and are therefore 
certainly filtering out the disk extended emission. 

Our best-fit results are presented in Table |2] and the cor- 
responding image is displayed in Fig. |2] We could not con- 
strain the kinetic temperature because the emission is opti- 
cally thin within the ring and we therefore assumed a con- 
stant temperature of 35 K, as found from the optically thick 



"CO J=2^1 emission by ' Guilloteau etall ( Il999l) . The de- 
rived inclination and p osition angle are in per fect agreement 



with those measured bv lGuilloteau et al.l (Il999h . Compared to 
the results of these authors, the ring width is narrower (50 AU 
instead of 80 AU) and consequently the surface density has 



to be higher (2.2 instead of 1.4 10 cm ) in order to repro- 
duce the observed flux density. The total mass derived from the 
model is 0.09 Mq. 

The ring width we measure remains only marginally re- 
solved at the angular resolution of our observations, while the 
ring edges have a very narrow width (Gaussians 1/e widths 
<10 AU), indicating a very sharp decrease in the density at the 
ring edges. Hence, the picture emerging is that of a very thin 
ring with a very pronounced overdensity compared to the inner 
and outer regions. 



4.2. The central source 



While Guilloteau et al. (Il999h noted the presence of emission 
within the GG Tau ring, the data presented in this paper clearly 
separate for the first time the ring from the central emission. 
The central source that we detect is unresolved at the resolu- 
tion of these observations (~50 AU). It is offset by 0.1" ± 0.1" 
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Ring physical parameters 



Inclination ( ) 


^7 -t- 1 


Position angle (°) 


7 ± 1 


Inner radius (AU) 


200 + 2 


Outer radius (AU) 


250 + 2 


Edge width (AU) 


< 10 


Temperature (K) 


[35] 


Temperature exponent 


[0] 


Surface density (cm""^) 


2.2+0.05x10^'* 


Surface density (g.cm^^) 


9.6+0.2 


Dust emissivity k (cm^'^/g) 


[0.012] 



Table 2. Fitted parameters of the GG Tau circumbinary ring. 
Error bars are 1 cr formal errors, not taking into account phase 
noise nor flux calibration uncertainty. Square brackets indicate 
fixed parameters. 



in R.A. and -0.1" + 0.1" in declination from the ring center 
position and has a flux density of 10 mJy, i.e. a brightness tem- 
perature of 1 .4 K in our synthetized beam. 



Using orbital parameters derived bv lBeust & Dutrevl(l2005b 
for the inner binary, we extrapolated the position of the stars at 
the time of the observations, i.e. beginning of 2008. These po- 
sitions (0.25" separation, -19° separation angle) are indicated 
in the top right panel of Fig. |2] The absolute astrometry of the 
stars not being known with enough precision, we forced the 
center of mass of the b inary (the two stars have similar masses, 
0.78 Mo and 0.68 Mq, ' White et al.l[l999h to coincide with the 
geometrical center of the ring (cross in Fig. |2). This hypoth- 
esis would not be coiTect if the ling were elliptical: the cen- 
ter of mass would then be at one of the ellipse focii. The data 
and models do however indicate the ring is circular at our an- 
gular resolution. With these assumptions, the inner continuum 
emission is clearly associated with the binary stars. It is thus 
tempting to identify this source with the thermal dust emission 
of one or two circumstellar disks around the individual stars 
(an alternative, although less probable origin would be free- 
free emission from an ionized jet). 

The continuum emission seems to peak at GG Tau Aa, al- 
though the formal fitted position indicates a marginal result 
(2.5 cr). Both stars exhibit optical emission lines that are typ- 
ical of classical T Tauri stars and trace accretion processes 
( White et al. 19991) . hence indicate that both sources should 



have an accretion disk. Our observations suggest that the mil- 
limeter continuum emission of the GG Tau Ab disk is weaker 
than that of the Aa disk, by at least a factor of 3, which points 
towards a less massive and/or colder disk. This is consistent 
with GGTauAb d isplaying weaker emission lines than Aa 
( White et al. 19991) . which is indicative of less intense accre- 



tion activity. 

A disk around GG Tau Aa would be truncated by the rel- 
ative motion of Ab to a size typically one third of the orbital 
distance, i.e. ~10 AU. Correcting for beam dilution, the bright- 
ness temperature can thus be ~30 K, which is a lower limit to 
the disk kinetic temperature Tkin- The coiTesponding mass, as- 
suming optically thin emission, would be 1 .5 lO"-' M© (scaling 
with 30 K/rki„). 



4.3. Streamer? 

The image of GGTau presented in this paper also reveals a 
continuous structure connecting the northern part of the cir- 
cumbinary disk with the central source. It is very suggestive of 
a streamer feeding matter onto the central individual disk(s), 
as suggested by a num ber of theoretical models of circumbi- 
nary disk evolution (e.g. lArtvmowicz & Lubowll996h . We note 
however that this structure is only marginally detected in our 
observation (contours are 2 cr in Fig.|2|i so that it is difficult to 
draw firm conclusions about its existence. 

Fig. |2] also shows several other low-level (one contour) ex- 
tensions of the ring. The reaUty of these structures remains to 
be confirmed, since they are close to the noise level and may 
also correspond to deconvolution artefacts in the imaging pro- 
cess. We note, however, that the position of the feature at the 
external edge of the south-east portion of the ring coincides 
with th e ' kink' seen in N I R sca ttered emission by Silber et alJ 
(l2000h 



IVIcCabe et al. ( 20021) . Similarly, the western-most 



extension is seen with several contours, and may trace a local 
enhancement of the column density, which could explain the 
west ern part of the ring bei ng brighter at lower frequency (3 .4 
Guilloteau et"aDll99 9'or 2.7 mm.'Putr ev et al.lll994l) . 



mm 



Finally, the apparent clumpy structure of the ring is not sig- 
nificant: all peaks correspond to a < 2cr difference with a uni- 
form ring, as for instance in the model shown in the right part 
of the figure. 

5. Summary 

We have presented a map of the 267 GHz continuum emis- 
sion in GG Tau A, at the unprecedented angular resolution of 
0.45 X 0.25", i.e. 65 x 35 AU at the distance of Taurus. These 
observations confirm the previous findings of Guilloteau et al. 
(1999) but allow us to place more accurate constrains on the 
properties of this system. Our main conclusions are: 

- The circumbinary dust emission is confined to a circular 
ring of width <50 AU, with very prononced sharp edges. 
These observations place very tight constraints on theoreti- 
cal model of the dynamical processes leadi ng to the forma 



tion of rings in circumbinary systems (e.g lBeust & Dutrev 



2005, 200 ^. GGTau is 1.5 Myrs old, but, as noted by 
Beust & Dutrev (I 061 observed ring may be a tran- 
sient feature that will undergo significant evolution in the 
next IVlyrs. 

- The data very clearly diffentiate the circumbinary disk from 
the emission located within the disk: an unresolved source 
is detected at the position of GG Tau Aa. This could be the 
first direct detection of of a circumstellar disk surrounding 
one of the stars of the central binary system. The mass of 
the circumprimary disk(s) is 1.5 lO"-' M©, about 80 times 
less than the mass of the circumbinary material. 

These observations confirm that GG Tau A is one of the best 
known and most interesting examples of a disk surrounding a 
young binary T Tauri system. It will be a piime target for future 
high angular studies with an instrument like ALJVIA. 
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